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1.  INTRODUCTION 

In  the  area  of  signal  processing,  a  model  that  is  often  used  involves 
expressing  the  observation  vector  as  the  sum  of  noise  vector  and  vector  of 
linear  combinations  of  the  (random)  signal  vector.  The  noise  vector  and 
signal  vectors  are  usually  assumed  to  be  distributed  independently  as  normal 
with  zero  mean  vectors.  When  the  noise  is  white,  the  problem  of  detection  of 
the  number  of  signals  transmitted  is  related  to  finding  the  multiplicity  of 
the  smallest  eigenvalue  of  the  covariance  matrix  of  the  observation  vector. 

So,  eigenvalue  methods  play  an  important  role  in  signal  processing.  These 
methods  play  a  dominant  role  in  the  area  of  multivariate  statistical  analysis. 

Some  workers  (e.g. ,  see  Kumaresan  and  Tufts(1980),  Li ggett( 1973) ,  Schmidt( 1979) , 
Tufts,  Kirsteins  and  Kurmaresan(1983) ,  Wax  and  Kailath(1984))  in  signal  processing 
have  used  the  eigenvalue  methods. 

Recently,  eigenvalue  methods  involving  information  theoretic  criteria  are 
used  by  Wax  and  Kail ath  (1984)  and  Zhao,  Krishnaiah  and  Bai  (1985)  for  deter¬ 
mination  of  the  number  of  signals  in  presence  of  white  noise.  The  object  of 
this  paper  is  to  detect  the  number  of  signals  present  in  presence  of  colored 
noise.  This  problem  is  equivalent  to  the  problem  of  studying  the  rank  of  r  when 
s2  =  r  +  ASj,  a  is  a  known  or  unknown  scalar,  ^  and  Zj  are  pxp  covariance  matrices 
and  r  is  nonnegative  definite  matrix  of  unknown  rank  q  <  p.  This  problem  arises 
in  other  areas  like  one-way  multivariate  components  of  variance  model  and  factor 
analysis.  Now,  let  n^  and  n2S2  be  distributed  independently  as  central  Wishart 
matrices  with  n1  and  n2  degrees  of  freedom  and  E(S.)  =  z  • , ( i =1 ,2 ) .  Rao  (1983) 
derived  the  likelihood  ratio  test  ( LRT)  statistic  for  the  rank  of  r  when  A  is 
unknown.  He  also  derived  a  modified  LRT  statistic  for  the  rank  of  r  when  A  is 
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known.  The  main  contribution  of  our  paper  is  to  propose  certain  information 
theoretic  criteria  for  detection  of  the  number  of  signals  and  establish  the 
property  of  strong  consistency.  The  paper  is  organized  as  follows. 

In  Section  2  of  the  paper,  we  discuss  the  model  considered  in  the  case  of 
colored  noise.  In  Section  3,  we  discuss  the  LRT  and  other  test  procedures  for 
testing  the  hypothesis  that  the  last  few  eigenvalues  of  are  eclua^  to  X  for 

the  cases  when  X  is  known  and  unknown.  We  propose  certain  information  theoretic 
criteria  in  Section  4  for  detection  of  the  number  of  signals  when  x  is  known  as 
well  as  when  it  is  unknown  and  establish  the  strong  consistency  of  these 
procedures.  Some  alternative  information  theoretic  criteria  are  also  mentioned. 

In  Section  5,  we  discuss  the  applications  of  the  above  results  to  determine 
the  rank  of  the  covariance  matrix  of  random  effects  vector  under  multivariate 
one-way  classification  model.  The  results  in  Sections  3-5  are  discussed  when 
the  first  q  eigenvalues  of  are  simple  (distinct).  When  these 

eigenvalues  have  multiplicities  (that  is,  they  are  equal  in  groups),  the  situation 
becomes  complicated.  In  this  case,  the  problem  invlolves  not  only  estimation 
of  q  but  also  the  multiplicities  of  the  first  q  eigenvalues.  This  problem  is 
investigated  in  Section  6. 
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2.  A  MODEL  IN  SIGNAL  PROCESSING 


In  the  area  of  signal  processing,  the  following  model  is  used: 


x ( t )  =  As(t)  +  xn(t)  (2.1) 

where  x(t):  pxl  is  the  observation  vector,  A  =  [A($j) ,.. .  ,A($  )] ,  s'(t)=(s1(t) 

•  ••>s  (t)),  s^t)  is  a  complex  waveform  which  is  referred  to  as  i-th  signal, 

A($^)  is  pxl  complex  vector  which  depends  upon  the  vector  <t>.  of  unknown 
parameters  associated  with  i-th  signal,  n(t)  is  a  complex  vector  associated 
with  the  noise  and  A  is  known  or  unknown  scalar.  We  assume  that  s(t)  and 
n(t)  are  distributed  independently  as  complex  multivariate  normal  with 
covariance  matrices  <1*  and  Ej  respectively,  E(s(t))=  0,  and  E(n(t))=0.  Also, 
A  denotes  the  complex  conjugate  of  A  and  A*  denotes  the  transpose  of  A.  The 
number  of  signals,  q,  transmitted  is  equal  to  the  rank  of  AfA*.  If  A,>. ..>\ 

1-  -  p 

1 

denote  the  eigenvalues  of  ,  q  is  given  by 


Ai  > 


>  A„>  A 


q+i 


(2.1) 


since  ^  ~  AtA  +  ae^.  When  E^  -  I,  the  problem  of  determination  of  the 

number  of  signals  was  considered  in  the  literature.  Wax  and  Kailath  (1985) 

used  Akaike's  IAIC  criterion  and  the  minimum  discription  length  (MDL)  criterion 

o 

due  to  Ri ssanen  and  Schwartz  when  o  is  unknown,  A=1  and  the  underlying  distribution 

is  complex  multivariate  normal.  Zhao,  Krishnaiah  and  Bai  (1985)  considered  an 

alternative  criterion  and  established  its  strong  consistency  for  the  cases 
2 

when  a  is  known  and  unknown  when  the  underlying  distribution  is  complex 
multivariate  normal.  They  have  also  considered  certain  cases  when  the  under¬ 
lying  distribution  is  not  necessarily  complex  multivariate  normal.  But, 
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it  is  not  realistic  to  assume  always  that  the  noise  is  white.  We  assume 
that  the  covariance  matrix  E^  of  n(t)  is  arbitrary  and  an  independent  estimate 
S1  of  Xj  is  available  from  a  different  data  set.  Also,  we  assume  that  n 
independent  observations  x(tj) . .  ,x(tn)  are  available  on  x(t).  In  this  case, 

n2 

we  can  estimate  X'2  with  S2  where  n2S  =  J  x(t.)x*(t.).  Since  S,  and  S9  are 

1  c_ 

distributed  independently  as  complex  Wishart  matrices  with  n^  and  n2  degrees 

of  freedom  respectively,  E(S1)  =  Zj  and  E(S2)  =  x2  =  AfA*  +  ^  ,  the  methods 

developed  in  this  paper  are  useful  in  finding  the  number  of  signals  transmitted. 

We  will  develop  the  methodology  for  finding  q  such  that  A,  >  ...  >  a 

1  q 

>  Aq+1  =  Ap  =  A  where  Aj  _>  . . .  _>  Ap  are  the  eigenvalues  of  E^j1,  n.^ 

and  n2S2  are  distributed  independently  as  real  central  Wishart  matrices  with 
nl  and  n2  degrees  of  freedom,  E(Sj )  =  l.(i=l,2),  and  e2  =  A^A1  +  AEj.when  A 
is  a  real  scalar.  Here  A:  pxq  is  a  real  unknown  matrix,  f:  qxq  is  a  real, 
positive  definite  matrix.  The  above  methodology  needs  only  trivial  modifica¬ 
tion  when  njSj  and  n2S2  are  complex  Wishart  matirces. 
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3.  TESTS  FOR  THE  EQUALITY  OF  THE  LAST  FEW  EIGENVALUES  OF  L'^"1 

Let  rijSj  and  n^S^  be  distributed  independently  as  central  Wishart 
matrices  with  n^  and  n^  degrees  of  freedom  respectively,  E(S^)  = 

e(S2)  =  ^2  ancl  z2  =  +  XEl*  T^e  ^°9  likelihood  function  L(e)  is  given 

by 

2L(o)  =  -  njloglSjl  -  n2log| e2 |  -  njtrz"^  -  n^rz^.  (3.1) 


Let  Hk:  Aj  >_  ...  _>  >  Ak+1  -  ...  -  a  =  1.  We  first  calculate  sup  L(e)  where 

6£0k 

°k  is  the  paramteric  space  under  H^.  Let  the  eigenvalues  of  S2S^  be  > 

...  >  6p.  With  probability  one  we  have  fij  >  s2  >  . . .  >  &  >0.  We  know  that 

there  exists  two  nonsingular  matrices  R  and  R  such  that 


*1  =  RR' ,  z2  =  RAR'  (3.2) 

Sj  =  RR' ,  S2  =  RAR 

where  A  =  diag(xi . .  ,Ap)  and  a  =  diag(61 ,. . . ,6  ).  Let  R_1R  =  V  .  Then 

2L(e)  =  -  (n1+n2)log|RRl  |  -  n£l  og(x  j . . .  xp)  +  L^V'.A),  (3.3) 

where 

LX(V^A)=  (n1+n2)log|V'V|  -  n^rV'V  -  n2tr(A-1VAV ') .  (3.4) 


First  we  fix  A  =  diagUj ,. . . ,x  )  and  compute  Supy.L^V'  ,A).  If  we  take  partial 
derivative  of  L  with  respect  to  V1 ,  we  obtain  the  following  optimizing  equations 

2(n1+n2)V"1  -  2n1Vl  -  2n2.AV,A"1  =  0 


i .  e. 
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anvv  +  enAV*A'1V=  Ip, 


(3.5) 
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where  I  is  the  pxp  identity  matrix,  and 


an  =  ni/n’  Bn  =  n2^n  ’  n  =  ni  +  n2' 


(3.6) 


From  (3.5)  it  follows  that  AV'A  1V  is  symmetric  and  hence  AV,A-1V'=V,a“1Va.  Since 

•i 

*1  >  6  2  >  •••  >  6P  >  °»  AV'A”  v  is  diagonal  so  that  by  (3.5)  V’V  is  diagonal. 
Thus  there  is  an  orthogonal  matrix  Q  and  a  diagonal  matrix  D  =  diag[d  ,...,d  ], 

dl  -  d2  -  *  •  •  -  dp  >  0  such  that  v  =  QD.  Since  A  and  av'A"1^/  are  diagonal, 
av  A  V  -  A  v'A  VA2  =  a2dq'A  ^DA-  so  that  Q'A  is  also  diagonal  and  the 
diagonal  elements  are  the  same  as  those  of  A"1.  Again  by  (3.5)  we  know  that 
the  diagonal  elements  are  arranged  according  to  the  increasing  order.  Hence 
Q'A  =  A  1  or  equivalently. 


a-1q  =  Qa"1 


Substituting  this  into  (3.5)  we  find 


(3.7) 


and 


Ip  =  anV'V  +  AVA^V  -  a^'V  +  3  V'A"aVA 

,..-1 


i  • 


anV'V  +  PnDQ,A~'LQDA  =  anV'V  +  $  DQ'QDA^A 


=  VV(anI+3nA'iA) 


-1,-1 


|V'V|  =  JanI  +  3nA_1A| 


=  n 


1 


1"1(Wil6i) 


(3.8) 

(3.9) 


Also,  we  have  by  (3.5) 


-  r^trV'V  -  n2tr(AV'A"1V)  =  -  (n^n^p. 

By  (3.4),  (3.9)  and  (3.10) 

SuP|Ll(V,,A)=  (n1+n2){-_p0g(an+6nA:16i)  -  p} 


(3.10) 


So , 


(3.11) 
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Sup  2L(o)  =  Sup  n{ - 1 og | RR 1 1  -  g  log( A, ...  A.  )  -  p 

0^0^  Aj>_.  .  k 

P  k  (3.12) 

-  14+1l09(VBn!1)  -  1411°9<Venxil61-),> 


where  0^  denotes  the  parametric  space  when  is  true.  Let  t  =  #{i  _<  p:  6.  >:  1} 
Also,  let  d  =  min{k,r},  and  set 


\  =  n .I  {-log(an+enAT16.)  -  pjogx.} 

1-1  (3.13) 

k  , 

*?  =  n  I  {-log(a  +3  A,  6.)  -  p  logA.}. 
c  i=d+l  n  n  1  1  n  i 


We  note  that  6^ 
for  i=l,. . .  ,d. 


>  >  . . .  >  5  .  >  1 ,  and  Sup  (j>,  can  be  reached  at  A.  =6 

Ai>...>A(.>l  1  1 

For  i=d+l,...,k,  6..  <  1  and  A^  >  1,  we  see  that  the  function 


VV  = 


log(an+Vil6i)  "  enlogXi 


(3.14) 


has  negative  derivative,  and  f. (a-)  is  decreasing  and  continuous. 

Sup  f • ( A . )  =  f.(l),  i=d+l ,. . .  ,k. 

A .  >1  1  1  1 

1 

From  the  above  discussion,  we  have 

d 


Sup  (L  =  -  ne  l  logs. , 
A1>...>A(J>1  r  i=l  1 


Thus 


(3.15) 


and 


Sup 

vr- 


*2  '  ~ 


k 

n  I  log(a  +6  6,). 
i=d+l  n  n  1 


(3.16) 
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From  (3. 12)  ,(3. 13)-(3. 15) ,  it  follows  that 

P 

Sup  2L(f>)  =  -  n{ log j  RR*  |  +  p  +  3  £  logfi.} 
e£©k  ni=l  1 


-  n 


P 

l 

i=l+inin(k,T) 


[1°9(an+Vp  -  yogs.] 


-  n^logi Sj|  -  n2log|S2|  -  np 


(3.17) 


n  £  ,  nogu  +3  y  -  ^nlo96^ 

i=l+min(k ,t)  n  n  1  n  1 


So,  the  LRT  statistic  for  testing  the  hypothesis  Hk  against  the  alternative 
that  the  rank  is  more  than  k  is  given  by 


Lk  = 


P 

n 


i=l+min(k,T) 


(a  5 

n  n  i 


)-n/26nBn/2 


(3.18) 


The  LRT  statistic  for  testing  against  (k<t)  is  given  by 


'kt 


min(t  ,t ) 

n 

i=l+mi n(k,x) 


/  „  ,-n/2  nfn/2 


n  n  i 


(3.19) 


if  k  <  t.  If  k  > 


then  L.  =  1. 


Rao  (1983)  considered  the  problem  of  testing  the  hypothesis  that  the 
rank  of  r  is  k  against  the  alternative  that  it  is  greater  than  k  when  = 

r  +  ae.  and  r  is  nonnegative  definite  for  the  cases  when  A  =  1  and  when  A  is 

unknown.  He  proposed  a  modified  LRT  procedure  and  the  LRT  procedure  for  testing 

the  hypothesis  on  the  rank  of  r  according  as  A  =  1  and  A  is  unknown. 

~k 

When  A  is  unknown,  let  denote  the  hypothesis  that 
A1  -  **•  -  Ak  >  Ak+1  =  =  Ap  =  A 
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for  k  -  0,1 . ,(p-l) .  Let  denote  the  model  for  which  is  true.  It  is 

known  (see  Rao(1983))  that  the  supremum  of  the  logarithm  of  the  likelihood 
* 

function  under  Hk  is  given  by 

-  -j-  1  og |  SJ  -  log|  S2|  -  ^ 

P 

+  *  i=Ll[anl09Ak0+3nl096j  “  l09(  Vko  +Vj  )]  (3.20) 


where  AkQ.  satisfies  the  equa'ti 


on 


P  6. 

p-k  =  S  - J- 


J=c+1  Vko+Vj 


or  equivalently. 


p-k  - 


ko 


j=k+l  a  A.  +B  6. 

n  ko  n  j 


(3.21) 


(3.22) 


As  pointed  out  in  Rao  (1983),  the  logarithm  of  the  likelihood  ratio  statistic 
is  given  by 


-  21og  Lk  =  log  ii  [((n?6.+n,A.  )/n)' 
i =k+l  6  1  1  K0 


1 


1  ko 


(3.23) 


which  is  distributed  as  chi-square  with  [(p-k) (p-k+l)-2]/2  degrees  of  freedom 
as  n^  and  n2  tend  to  infinity. 

We  will  propose  the  following  alternative  procedures  for  testing  the 
hypothesis  Hk  against  the  alternative  that  Ak+^  >  1.  We  accept  or  reject  Hk 
against  Ak+^  >  1  according  as 
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g(£k+l’“-  ’V  >  Cc 


(3.24) 


where 


P[gUk+i’*--’'U  Ji  “  (!“<*)• 


V  -  VV  -  ^~a)‘  (3.25) 

For  example,  g^,...,^)  may  be  £k+1  or  £k+1+  ...  +  £p.  The  exact 
distributions  of  the  above  statistics  are  complicated.  Also,  they  involve 
nuisance  parameters  unless  k  =  0.  But,  the  joint  asymptotic  distribution  of 
is  9i ven  *n  a  companion  paper  (in  preparation)  by  the  authors 
for  the  real  and  complex  cases.  A  review  of  other  asymptotic  results  was 
given  in  Murihead  (1978).  We  can  use  the  above  result  to  obtain  asymptotic 
distributions  of  statistics  like  l.  and  L  +...+£  . 

K  “  1  K+i  p 

We  will  now  consider  the  case  when  X  is  unknown.  Let  H..:  A.  =  x..  Then 
Hk*can  be  deceased  as  ft  iV  ft  ^  >w  and  .(^^y  Motivated  by  the 


above  decompositions,  we  propose  the  following  procedures.  We  accept  H 


P-1 

against  [x.>xj  if 
i=k+l  1  P 


<W  -  C«1 


for  i  k+1 , . . .  ,p-l  and  reject  it  otherwise  where 


(3.26) 


p[(^k+l^p)  1  cal^Hk^  ~  (1-“)‘ 

-k 

If  Hk  is  rejected,  we  accept  or  reject  the  subhypothesis  H* 


(3.27) 


accordinq  as 


IJl  >  c 


1  P  a  1 

The  hypothesis  H'  when  tested  against  (j  t>,  >  A,+1 )  is  accepted  if 

i=k+l  1  1+1 


for  i =k+2 , . . . ,p-l. 

* 


(3.28) 


for  i  -  k+1 ,2 ,...,( p-1 )  and  rejected  otherwise  where 


P[(£j/£.j+i)  £  ca2>  i  k+1 ,2,...  ,p-l|H.  ]  -  (1-a). 


(3.30) 


Similarly,  the  hypothesis  H.  when  tested  against  [a -  >  X . , .  ]  is  accepted  if 

i<j=k+l  1  1+1 

and  only  if 


for  i  <  j  and  i  =  k+1,..., p-1  where 

P[(£k+l/£p)  -  Ca3lHk]  =  (l'a)'  (3-31) 

We  can  also  use  '^|<+]/(£|<+]+*  •  -+^p)  as  a  test  statistic.  Exact  computations  of 
the  probability  integrals  associated  with  the  above  procedures  are  complicated 
and  involve  nuisance  parameters  except  when  k  =  0.  In  this  particular  case, 
percentage  points  are  available  for  a  few  special  cases  (see  Krishnaiah  (1980)). 
But,  approximations  to  the  critical  values  c  . ,  c  0  and  c  n  can  be  derived 

al  a<i  ao 

for  large  samples. 
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4.  DETECTION  OF  NUMBER  OF  SIGNALS  USING 
INFORMATION  THFORETIC  CRITERIA 


In  the  preceding  section,  we  discussed  various  procedures  for  testing 
the  hypotheses  on  the  number  of  signals.  We  will  now  discuss  procedures  for 
detection  of  the  number  of  signals  by  using  information  theoretic  criteria. 
When  A  =  1,  the  likelihood  ratio  test  statistic  Lk  for  Hk  is  given  by 

G(k)  =  logLk 

I (k »Cn)  =  -  logLk  +  Cnv(k,p)  (4.1) 


(3.18).  Now,  let 
and 


where  v(k,p)  =  ^-k( 2p- k+1 )  denotes  the  number  of  free  parameters  and  C^ 
satisfies  the  following  condition: 


(i)  lim(C  /n)  =  0 

n~*oo 


(4.2) 


(ii)  lim(Cn/loqlogn)  =  (4.3) 

n-*» 

A 

Then,  according  to  FDC  criterion  we  find  q  which  satisfies 

I(5,Cn)  =  min{I(0,Cn),...,I(p-l,Cn)}  (4.4) 

and  use  q  as  an  estimate  of  q  which  is  the  number  of  signals  present  in  the 
true  model  M  ,  The  strong  consistency  of  q  is  proved  below: 

THEOREM  4.1.  If  ~ W  (n.  ,z.) ,  i=l,2,  n  =  nx+n2  +  »  and  an£[a ,b]  d  (0,1) 

A 

with  a,b  being  constants,  then  q  is  a  strongly  consistent  estimate  of  q  under 

the  model  M  . 

H 

PROOF.  It  is  known  (see  Zhao,  Krishnaiah  and  Bai  (1985))  that 


I*,  -  ql  - 


a.  s. 


(4.5) 


13 


Using  Taylor's  expansion,  we  get 

log(an+fW  ‘  Bnlog6i  =1  Vn(V1)2(1+o(1))  a.s.  (4.6) 

for  i  >  q.  Here  we  used  the  fact  that 

lim  6.  =  1  a.s.  for  i  >  q.  (4.7) 

n-H*. 

With  probability  one  for  large  n,  we  have  min(q,T)  =  q  so  that,  for  large  n, 

|G  (q)  -  G  ( k) |  =  O(loglogn)  =  o(Cn)  (4.8) 

and 

(Kq,Cn)  -  I(k,Cn))/Cn  =  "  (G(q)  -  G(k)/Cn  -  (k-q)  (2p-k-q+l)  (4.9) 

-»■  -  ( k-q)  (2  p-k-a+1) 

when  k  >  q.  Thus  with  probability  one,  for  large  n, 

I(q,Cn)  <  I(k,Cn)  ,  if  k  >  q.  (4.10) 

Since,  with  probability  one,  min(q,T)  =  q,  for  large  n,  we  have 

q 

G  (q)  -  G  (k)  =  n  £  [log(a  +3  6.)  -  BjogS.]  (4.11) 

i=k+l  n  n  1  n  1 

for  k  <  q.  Note  that  for  i  _<  q ,  lim  6.  =  A.  >1,  we  see  that  there  exists 

n-*»  1 

a  constant  y  >  1  such  that  with  probability  one 

<5.  >  y  for  i  =  k+1 . q, 

for  large  n.  By  the  monotonicity  of 

f(fi)  =  log(an+Bn<5)  -  3nlog6, 


(4.12) 
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P ( CD)  =  P[I(q,Cn)  -  I(k,Cn)  <  0;  k=0,l ,. . . ,(p-l) ;  Mq|Hq].  (4.20) 

It  would  be  of  interest  to  compare  numerically  the  probability  of  correct 
detection  of  the  criteria  I(q,CN),  AIC(q)  and  MDL(q). 

We  now  discuss  the  problem  of  detection  of  the  number  of  signals  when  A 
is  unknown.  In  this  case,  the  logarithm  of  the  LRT  statistic  for  H*  is  given 
by  (3.23).  Now,  let 


G  (k)  ~  logl^ 


(4.21) 


and  assume  that  Cn  satisfies  the  conditions  (4.2)  and 

(4.3).  Then  an  estiamte  of  q,  the  true  number  of  signals,  is  given  by  q  where 


q  =  max{k:  l<k<p-l,  G*(k)  -  G*(k-1)  >  Cn}  (4.22) 

and  maxtfi  =  0  for  convenience. 

★  * 

Let  Mk  denote  the  model  under  which  Hk  is  true.  We  now  prove  the  strong 
consistency  of  the  above  method. 


THEOREM  4.2.  If  S.  ~ 

-  i 

a,b  being  constants,  then, 
estimate  of  q. 


^p(ni,Ei)’  1  ~  1»2»  n  ->■  and  ane[a,b]  e(0,l)  with 

* 

under  the  true  model  Mq ,  q  is  a  strongly  consistent 


PROOF.  As  pointed  out  earlier 


6.  -  A. 


=  O(^loglogn) 


a.s 


(4.23) 


for  j  =  1 ,. . .  ,p. 

Note  that  bei ng  6^  >  62  >  . . .  >  with  probability  one,  we  can  see  from 

(3.21)  or  (3.22)  that  Ako  >  Ak+1  Q  for  0  <  k  <  p-1. 

* 

We  assume  that  M  is  the  true  model,  and  k  q.  As  mentioned  above, 

M 
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l6j  -  1|  =  0(\j^logTogn)  a.s. 
for  j  >  q.  Assume  I «  •  -  1 1  >  ep  for  j  >  k :+  1.  Then 


2d< 


P 


6  . 


_  /  — - — j - =  p-k  <  — --i-P.-k) - 

“nW1"6^  +Bn  J=k+1  anAkO+pn6j  anAk0^1+en^"1+Bn 

A 

and  1  -  <_  Ak0  _<  1  +  en>  Thus  we  have 


Ak0  "  =  °(\  hr'oglogn)  a.s.,  k  >  q. 


>  q. 


Using  Taylor's  expansion,  we  see  that  for  k  _>  q, 

P 

Vll[“nl°9  k0  +  6nl°94j  "  109<“n'kO+lin5j)] 

P 

7  j=Ll[‘a"<X"k0'1)2  '  6"(6i'1)2  +  ‘“n^ko-1)  +  e„(«ri)): 


a.s.  n 


=  0(  loglogn) 


a.s. 


If  k  >  q ,  then  from  lim  C^/  loglogn  =  °°,  we  get 

n-*» 


L*(k)  -  L*(k-1)  =  o(Cn)  a.s. 

Thus,  with  probability  one,  we  have  for  large  n, 

L*(k)  -L*(k-1)  <  Cp  for  k  >  q  , 

A 

which  implies  q  =  0  if  q  =  0.  Now  we  assume  that  1  <_  k  <  q.  We  have 

G  (k)  - G*(k-l)  =  n(g(xk0)  -  g(x“k.i,0)  -  %Wk.lj0  -  V°9sk 
+  lo9<Vk-l,0+6n{k»- 


(4.24) 


(4.25) 


](l+o(l)) 


(4.26) 


(4.27) 


(4.28) 
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where 


g(x)  -  a  (p-k)logx  -  l  log(a.nx+B  <s.) 

i-Cxi  •  •  n  1 


j=k+l 

It  is  easily  seen  that  for  X£(^|<q»^|<  j  q]  , 

9'(X)^"[:I  W7-(p-k)]< 

So  we  have 

9(\0)  -  g(\_1}0)  >  0. 

Let  pn  be  the  solution  of  the  following  equation: 

P  A . 

p  -  k  +  1  =  l  — J - 


j=k  Vn+3nY 


(4.29) 


(4.30) 


(4.31) 


From 


p  -  k  +  1 


6  . 


J  k  anAk-l,0+3n6j 


and  lim  6.  -  x.  a.s.,  it  follows  that 

n->oo  J  j 


^Vl.Cf-n) 


li 

n-*» 


-  0  a.s . 


(4.32) 


Thus  we  have 


n{log(Vk-l)0+Vk)  "  anlogAk-l,0  "  3nlog<sk} 


a.s. 


n{lo9^Vn+BnV  "  anlogyn  "  3n1ogAk}{1  +  °(1)}  (4.33) 


Si  nee 


P 

A 


A  . 

JL 


“nWj 


=  1'  + 


A  . 
J 


j=K+l  Vk+3nAj 


<  1  +  P  -  k. 
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it  follows  Pp  <  yQ  <  for  some  constant  yQ. 

Hence 

l09<  Vn+Vk>  -  “n,09un  ‘  en,09\  i 

(4.34) 

log(anyo+BnAk)  "  an1ogiJo  "  en1ogAk  =  h(ct)* 

Note  that  h(a)  is  continuous  for  ae[a,b]  where  $  =  1-a.  So  there  exists  a 

a0c[a,b]  such  that 

h(an)  >  h(ao)  =  1og(aoyo+B.Ak)  -  ctologyo  -  30logXk  >  0  (4.35) 

By  (4.33)  -  (4.35),  we  see  that  with  probability  one  for  n  large, 

n{log(Vk-l,0HVk)  "  anl!ogAk-l,0  ~  3nlog6k}  (4.36) 

>  \  nh(aQ). 

From  (4.28) , ( 4. 30)  , ( 4. 36)  and  C^/n  ->  0,  it  follows  that  with  probability  one 
for  n  large, 

L(k)  -  L(k-l)  >  MCn,  k  =  1,2,. ...q  (4.37) 

for  any  fixed  M  >  0. 

By  (4.27)  and  (4.37)  we  see  that  with  probability  one  for  n  large, 

A 

q =  q>  (4.38) 


and  the  theorem  is  proved. 
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5.  MULTIVARIATE  ONE-WAY  RANDOM  EFFECTS  MODEL 

In  this  section,  we  discuss  the  relationship  between  drawinq  inference  on 
the  rank  of  the  covariance  matrix  of  column  effects  in  one-way  multivariate 
random  effects  model  and  the  problem  of  finding  the  number  of  signals  discussed 
in  the  preceding  sections.  The  one-way  multivariate  random  effects  model  is 
given  by 


for  i  -  l,2,...,k  j  =  1,2,...  ,m  where  y  is  the  general  mean  vector,  cu  is  the 

vector  of  random  effects  of  i-th  column  and  x. .  denotes  the  j-th  observation  on 

1  J 

i-th  column  and  £ .  .  is  distributed  as  multivariate  normal  with  mean  vector  0 
~  '  J 

and  covariance  matrix  Ej.  Also,  is  distributed  independent  of  ..  as 
multivariate  normal  with  mean  vector  0  and  covariance  matrix  ip.  The  covariance 
matrix  of  x.  .  is  where  E^  =  4)  +  E^.  It  is  of  interest  to  test  the  rank  of 
ip  is  r.  If  the  rank  of  ip  is  r,  then  there  exists  a  full  rank  matrix  B:(p-r)xp 
such  that  Bip  =  0.  Testing  the  hypothesis  that  the  rank  is  zero  is  equivalent 
to  testing  the  hypothesis  of  no  column  effects.  If  ip  is  not  of  full  rank,  then 
we  can  take  advantage  of  this  knowledge  in  estimating  ip.  Anderson  (1984,1985) 
and  Schott  and  Saw  (1984)  have  independently  derived  the  LRT  statistic  for 
testing  the  tiypothesis  on  the  rank  of  ip.  Now  let  S^  and  Sw  respectively  denote 
the  between  groups  and  within  group  sums  of  squares  and  cross  products  matrices 
respectively.  Then  S^  and  S are  distributed  independently  as  central  Wishart 
matrices  with  (k-1)  and  (km-k)  degrees  of  freedom  respectively,  E^)  =  (k-1) 
(Ej+mip)  and  E(SW)=  (km-kjE^ 
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60  DETECTION  OF  NUMBER  OF  SIGNALS  WHEN 
EIGENVALUES  OF  HAVE  MULTIPLICITIES 


In  Sections  3  and  4,  we  discussed  the  problem  of  detection  of  the  number 
of  signals  under  the  assumption  that  the  nonzero  eigenvalues  of  the  matrix  AyA' 
are  distinct.  But  situations  arise  when  it  is  unrealistic  to  make  the  above 
assumption.  In  this  section,  we  consider  the  problem  of  detecting  the  number 
of  signals  and  finding  the  multiplicities  of  the  eigenvalues  of  AfA'.  We  will 
first  discuss  the  problems  of  finding  the  rank  of  r  when  the  underlying  distri¬ 
bution  is  real  multivariate  normal. 

p 

For  the  interval  [o  ,p] ,  there  exists  £  ^  ^  different  integer 

partitions  such  as  0=kQ  <  ^  <  ...  <  k£  =  p ,  l  =  1,2,... ,p.  We  denote  the 
set  of  all  such  partitions  with  K.  Let 


Hk  k  :  Ak  +1  =  ...  =  A.  =  c. ;  i  =  1,2  ,...,£ 
k1***k£  Ki-1  ki  1 


(6.1) 


where  Cj  >  . . .  >  c^  are  unknown  constants  and  x1  >  a2  >  . . .  >  xp  are  eigenvalues 
of  We  denote  the  corresponding  parametric  space  and  model  with 


°k1...k„  and  Mk  ...k  resPectively.  We  are  interested  in  selecting 
★ 

model  M 


the  correct 


5*  •  • 


using  information  theoretic  criteria  when  we  do  not  have  any 


knowledge  of  qp...,qr.  When  H*  k  is  true,  the  log-likelihood  function 

*  1  *  ‘  *  Z 


L  (e)  is  given  by 


-k 

2L  (e)  =  -  rijlogl Ej |  -  n2log|z2|  -  r^trz^Sj  -  n^rz”^.  (6.2) 


At  first  we  calculate  Sup 


0e0k  ■  k 

9  •  •  • 


2L  (e). 


Denote  the  eigenvalues  of  S2S^  by  <5^  >_  .„.  _>  5p.  There  exists  two  non¬ 
singular  matrices  R  and  R  such  that 
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Ej_  =  RR' ,  l2  =  RAR' , 

A  /V  A  A 

Sj  =  RR ' ,  S2  =  RAR' , 

where  A  -  diag(Aj,. . .  ,Ap)  and  a  =  diag(<5^,. . .  ,s  ) .  Without  loss  of  generality, 
we  assume  that  >  . . .  >  sp  >  0.  Now,  let  R_1R  =  V.  Then 

2L  (e)  =  -  (n1+n2)log|  RR'  |  -  n2log(x1...xp)  +  L^V  ,a),  (6.3) 

where 

L1(V  ,A )  =  (n1+n2)log|V'V|  -  n^rV'V  -  n2tr( V 'A_1Va) .  (6.4) 

First  we  fix  A  =  diag(A^ ,. . . ,xp)  where  (6.1)  is  satisfied.  For  given  A, 
we  can  now  calculate  Supv  L1(V  ,A).  To  accomplish  this,  we  take  derivative  of 
with  respect  to  matrix  V1  and  obtain  the  following  optimizing  equations: 

2(n1+n2)V-1  -  Zn^'  -  2n2AV'A~1  =  0 

i  •  6  • , 

«nV'V  +  BnaVA-1V  -  Ip,  (6.5) 

where 

an  =  nl/n’  Bn  =  n2/n’  n  =  ni  +  n2*  (6.6) 


Using  the  same  argument  as  used  to  prove  (3.8),  we  find 

rV(“nIp+6nA'l4)  =  V 


and 


=  a 


-1  i -1  P 

„W  A|  -  H 


■  .  a  A .  +  3  <5 . 

1=1  n  i  pn  i 


(6.7) 


(6.8) 
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Also,  we  have  by  (6.5) 


-  rijtrV'V  -  n2tr(AV'A_1V)  =  -  (n^Jp. 


(6.9) 


Now,  using  (6.3),  (6.4),  (6J8)  and  (6.9),  we  get 
Sup 


0t0k  k 

>*  •  *  5^, 


2L  (e)=  Sup 

v*>ci 

l 


l 


n{-  1 og | RR 1 |  -  p  +  a  J  Ak.logc.  (6.10) 

i=l  1  1 


l  .  1  log(«  c.  +36.)} 

i=l  j«K,  n  1  n  J 


where  k..  -  {k^j+l,  k._i+2 ,. . .  ,ki }.  Now,  let  c..  's  be  chosen  such  that 


JeKi  anC. 


Aki  ,  i 


1  l 


(6.11) 


Noticing  that  c..  _>  <5^  >  6^  +1  _>  c^j  ,  then  we  get 


Sup 
0€0  ' 


1  1 


kp. . .  ,k^ 


2L(e)  -  n{-  anlog|S1|  -  3nl og | S2 |  -  p 


(6.12) 


A  .  I  (a  log  C,  +6  logfi.-log(a  C.  +36.))}. 
1=1  jer.  i  n  j  fii  n  j 


Also,  we  have 


“VIP*  =  ^  “n'031'5!1  '  6n’“9|S2|  -  pj. 


Sup 

T. 


(6.13) 
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Now,  let 


ic 

L  (kj , . . . ,k^)  =  Sup 


e^e 


{ 2 Lt  e ) }  -  Sup  (2 Ct  e )  > . 


kj . ,k^ 


z1>o,e2>o 


(6.14) 


Then,  we  get 


L  (klf...  ,k£)  =  n  l  l  (a  logc. +3  log6.-log(a  c.  +Bn6.))  (6.15) 

1  =1  jtK  .  J  11  1  n  J 


Now,  let 


6  (ki». . . , k^ )  -  L  (kj 


(6.16) 


where  satisfies  the  following  conditions: 


(i)  lim( C /n)  =  0 

n->oo 

(ii)  lim( C  /loglogn)  =  »  (6.17) 

it**5 

Then,  we  estimate  (r,q1 ,. . .  ,qr)  with  (r ^  ,. . .  ,qp  where 


G  (q1»...»qp)  =  max  G*(k, ,k„)  (6.18) 

(k15...,k  f)£K  1  1 

whereic  =  { (kj ,. . .  ,k£) ,  kj  <  k2  <  ...  <  k£  =  p,  l  =  1,2,...  ,p). 

We  now  prove  the  strong  consistency  of  the  above  procedure. 

THEOREM  6. 1 .  Let  n^S^  anc*  n2^2  distributed  independently  as  central 
Wishart  matrices  with  n^  and  n2  degrees  of  freedom  respectively.  Also,  let 
E(Sj )  =  ^  (i  =  1,2).  Then  ( r,q^ .  ,q~)  defined  by  (6.16)  is  a  strongly  consis¬ 
tent  estimate  of  (r ,q^ ,. . .  ,q^)  when  n  ->  «>,  [a,b]c(0,l)  with  a  and  b  geing 

constants  and  the  true  model  M* 


q 


r 


•q 


r 
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FTOOF.  Suppose  that  M*  is  the  true  model  and  X,  >  >  x  are 

Hj* •  -Hr  1  —  "  *  •  —  p  c 

the  eigenvalues  of  z^1.  By  the  law  of  the  iterated  logarithm,  we  have 
s,  -  q  =0(Npl2M)  a.s.  as  n  - 
for  i  =  1,2.  Thus, 


-  hVi*  =  0(\|M22!i)  a.s. 
By  Lemma  3.2  in  (Zhao,  Krishnaiah  and  Bai  (1985)), 


as  n  ->■  co 


I  _  xi  I  “  0(\--^°9n)  a.s.  as  n  oo 

fori  =l,...,p.  Suppose  that  (kj,...,^)  3  Then 


(6.19) 


Sup 

0tOJ 


2Ll9)  _>  Sup 


k -i  9...  ,k 


l 


060* 


2Ll(0), 


1  **  ”  ,Jm 


so,  we  have 


L  >_  L  (jj .  ,jm) .  (6.20) 

Now  we  suppose  that  (k1#...,k£)  a  (q1#...#qr).  Then  l  >  r.  Write 
ci  =  {Vl+1’  Vl+2’***’V’  Aq.  =  a.  and  v.  =  (6.-Xi)/x..  Also,  put  J  = 

A 

-2 _ _2  for  i  =1 ,2, , . .  ,r.  Assume  (1-eJ  c.  <  «.  <  (l+ejc.  for  all  jeK..  Then 

i 


ai 


1  l 


6  . 
J 


‘rfi  ci1(1-£n)'1+en  **1  %  S'  +Vi 


ai 


(6.21) 


Si  nee 


(6.22) 


lci  '  cil  -  en  ci  for  1  " 
Using  (6.19)  and  (6.22),  we  have 


A. 


(6.23) 


as  n  -*  °°  for  i  -  l,...,r.  By  (6.15)  and  Taylor's  expansion,  we  get 

0  <  -  Llqx,...,qr)  £  -  H[q1 , . . .  ,qr)  (6.24) 

r 

=  n  I  l  [log(l+a  IJ.+B  V.)  -  a  log(l+y.)  -  3  log(l+v • )] 
i=l  je  k.  J  i  m  j 

=  f  X  X  K^Vj  +  (Vl+Bnvi)2](1+o(1))  «.s. 

1-1  JC!\j 

Now  using  (6.19)  and  (6.23),  we  have 

0  1  -  L*(q15...  ,qr)  =  O(loglogn)  a.s.  (6.25) 

as  n  ->  ».  By  (6. 16) , (6. 17) ,  (6.24)  and  (6.25),  with  probability  one  for 
large  n, 

G*(q2,,"*,9y.)  -  G  (kj,. ,  k^)  -  L  (q1 ,. . .  ,q^,)  -  L  ( k^  , . . . ,  k^) 

+  (l-r) Cn  >  0.  (6.26) 

Finally  we  suppose  that  (kp.  ..,k^)  is  a  partition  of  [o,p]  such  that  there 
exists  at  least  one  qt  satisfying  k-_x  <  qt  <  k.  for  some  i.  Define  a  new 
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partition  (1,2,...  ,qt~l ,qt+l ,  qt+2,...,p)  =  (jj ,. . .  ,jp_1) .  By  the 

and  (6.20) ,  we  have 

L  (^1  ’•  •  •  ’Oy.)  -  L-tkj  ,. . .  ,k^)  _>  L?tq !  »• . .  ,qr)  -  L*(  jj ,. . .  ,jp_1) 
Now,  let  Nt  =  {qt,  qt+1>.  It  is  easy  to  see  that 

*  #  .  .  a  A 

-  L  .)  =  n  i  [log(a  A+3  6.)  -  a  logA  +  6  log<5.1 

F  jc-N,  nnj  n  n  3  jJ 


where 


l 

jc-Nt 


ot„A+g  6 . 

n  pn  j 


=  2. 


Define  yp  such  that 


A. 

1 


anWj 


=  2. 


By  lim  <5.  =  A  .  a0s.  for  jfeN.  ,  we  have 

rn-oo  J  j  z 


lim(A-p)  =  o  a.s. 

n-^o  " 


By  (6.31),ct  >  ct+1,  and  the  condition  an6[a,b] c(p#l) ,  there  exists 

y  such  that 

0 


ct  > 


n  * 


n  =  1,2,... 


By  (6.29),  (6.31)  and  (6.32),  we  get 


p-T 


a.s. 


jL  [1°9(‘WVj)  -  “nl09«n 

3nlogAj]  +  o(n) 


-  n[log(Vn+Bnct  }  “  an1o%  -  3nlogct]  +  o(n) 


fact 


(6.27) 


(6.28) 


(6.29) 


(6.30) 


(6.31) 

a  constant 

(6.32) 
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-  n[1°9(Vo+Pn  ct  }  “  an1ogyo  '  Bn1og  ct  ]  +  o(n)  (6.33) 

The  function 

h(a)  =  1og(apo+6ct)  -  anlogpo  -  Bjog  ct  ,  (e=l-o)  (6.34) 


is  positive  and  continuous  forae[a,b].  So  there  exists  a  constant  aQ£[a,b]  such 
that 

h(an)  >  h(a0)  =  ■>og(aouo+3o  ct  )  -  aologu0  -  B0logct  >  0.  (6.35) 

By  applying  (6.33)  -  (6.35),  we  get,  with  probability  one 

*  ^ 

>•  •  •  2'nh(aQ)  (6.36) 

for  large  n.  From  (6.25)  ,(6.27)  ,(6.36)  and  lim  C  /n  =  0,  we  know  that, 

n-*» 

with  probability  one  for  large  n, 

★ 

G  »•  •  •  >qr)  -  . .  ,k^)  =  L*(q1 ,. . .  ,qr)  -  L*(k1 ,. . .  ,k^)  +  o(n) 

>  \  nh(aQ)  +  o(n)  >  0.  (6.37) 


From  (6.26)  and  (6.37)  it  follows  that,  with  probability  one  for  large  n, 

A  /\ 

(r,q15...,q-)  =  (  r^ , . . .  ,qr) ,  (6.38) 

which  completes  the  proof  of  Theorem  6.1. 

★ 

REMARK  6.1.  L  (k^,...,k^)  can  be  regarded  as  a  general  test  statistic  (not  necessarily 

LRT  statistic)  for  testing  the  hypothesis  H*  .  .  Also,  letruS-  and  n0S0  be 
nj  Kr**K£  11  22 

distributed  as  J  X..X.  and  £  Y.Y respectively  where  X19...,X  and  Y, ,. . .  ,Y 
i=l~  i=l~  -1  ~nj  ~1  ~np 
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are  subject  to  the  following  conditions: 

(a)  X15...,Xn^  are  i.i.d.  and  YJ..-..Y  are  i.i.d.  such  that  E(X1)=E(Y1)=0 

(b)  (X, . X  )  an,d  (Y1 .  ,Y  )  are  independent  or  not 


(c)  E(XjXj)  =  Zj  and  E ( YjY^ )  =  e£ 


are  positive  definite 


(d)  ECXjXj)2  <  co  and  E(Y^Y1)2  <  » 

A  /\  A 

Then,  (  r ^ . qp  is  still  a  strongly  consistent  estimate  of  (r.q^.. .  ,q  ). 

REMARK  6.2.  If  =  c^  is  known,  we  can  assume  c^  =  1.  In  this  case,  we 


redenote  H. 

Kr 


»k« 


and  G  (k 


1 


.y  by 


the  corresponding  notations  without  stars.  Following  the  same  lines  as  in  the 
case  of  proof  from  (6.1)  to  (6.10),  we  find  that 

n  n 

sup  2L  (e)  =  Sup  n{-  1 o  g [ RR|  +  a  V  logx  -  J  log(a  A.+BS, )p} 

0(kr..k£)  c1>...>c£_1>l  n  i=l  1  i=l  n  1  n  i 

l-l  l-l 

~  Sup  n{-  log | RR|  -p  +  a  ^Ak.logc.  -  l  l  log(a  c.+e  6. ) 

cl>. . . >c£_i>1  i=l  1  1  i=l  jfc«.  ni  ni 


-  .1  log(an+3n6i)} 


(6.39) 


Define  c^  as  those  defined  in  (6.11)  and  define 

t  =  maxfi  <_l-l,  c.  >  1} 


(6.40) 


If  i  <  t,  then 


Sup  (a  Ak.logc.-  I  log(a  c.+Mj)  = 
c.j>l  jcKi  n  1  n  1 

=  anAkilogci  -  .I  log(anc.+B  6i) 

J<=K.j 


(6.41) 
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and  the  superium  can  be  reached  at  c^  =  c..  If  i  >  T,  i.e.,  c^ 
we  have 

gf-I^Ak.logc,-  l  l09(Viten6i>) 


JtK.j 


-  qn&ki  v _ °n  %/'  v  ci  \ 

Ci  Vl+Vl  "  Vi+Vi ) 


U  / 


ci 


=  0, 


j6Ki  anci+Bn6v 


Hence 

SuM^AkJogc.j-  ^  pog(anci+Bn6i))  =  -  .  I  lo9(VVi 
and  the  superium  is  reached  at  =1.  Noting  that  Cj  >  ...  > 

i(kr..k,)2L*(9)  =  n<_  “"l09|Si1 '  e"l09|s2| ' p 

T  £ 

+  X  A  (anlogci+Bnlog6j'1og(Vi+en5i);i  “  .  I  1og(ctn+Bn 

1  JfcK-j  i  -T  +  l 

Hence  we  get 

T 

L(kr..k£)  =  n  l  l  [a  logc,+3  logs.  -  log(a  c.+b  S.) 

i  =  l  je  k .  '  mini 

l 

-  n  t  log(a  +BnS.) 
i  =t+1  n  n  1 

Now,  let 

G(kj,...,k^)  =  L(kj,...,k^)  -  £cn 

where  cn  satisfies  (6.17).  Then  we  estimate  (r,q^,...,q  )  with 
where 


<  1 ,  for  c^  >  1 


). 

c  >  1 ,  we  obtai n 

T 


!-)>  (6.42) 


(6.43) 


(6.44) 

A  A  A 

(r,q15...  ,qp. 
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A  A 

G(q-,  ,q;)  =  max  G(k- ,k») ,  (6.45) 

1  JL 

where  K  is  the  set  defined  below  (6.18).  Similarily  as  Theorem  6.1,  we  have 
THEOREM  6.2.  Let  r^S,,,  be  distributed  independently  as  central  Wishart 
matrices  with  n^  and  n^  degrees  of  freedom  respectively.  Also,  let  ES.  =  x.(i=l,2). 

A  A  A 

Then  (r,q^ ,. .. ,qp)  defined  by  (6.45)  is  a  strongly  consistent  estimate  of 

(r,q1 ,. . .  ,qr) ,  when  n  -*•  a^efa  ,b]  c (0,1)  with  a,b  being  constants  and  the  true 
model  M 

q  1  •  •  •  q  r 

Since,  with  probability  one,  when  large  enough,  x  =  r  -  1  under  the  hypothesis 

H  ,  we  find  that  Theorem  6.2  can  be  proved  by  the  same  argument  as  used 

U1‘ '  -r 

in  the  proof  of  Theorem  6.1.  So  that  we  omit  it  here. 

REMARK  6.3.  Remark  6.1  concerning  Theorem  6.1  is  also  available  to  Theorem  6.2. 
REMARK  6.4.  If  n^S^  and  n2S2  are  distributed  independently  as  central  complex 
Wishart  matrices.  Then  the  log-likelihood  function  L(e)  is  given  by,  up  to  an 
adding  constant, 

L(o)  =  -  n11og|  z1|  -  n2log|>:2|  -  n1tn:"1S1  -  n^rz"^. 

In  the  arguments  in  Sections  3-6,  we  only  need  change  the  following  notations 

=  RR*,  z2  =  RAR* 

S^  =  RR  S2  =  RAR  in  (3.2) 

-k 

where  A  denotes  the  transpose  of  the  conjugate  of  the  matrix  A. 

4(V',a)  =  (n1+n2)log|V*V|  -  njtrV*V  -  n^rA^VAV*  in  (3.4) 

and  Q  being  a  Hermite  matrix  instead  of  an  orthogonal  matrix  and  rewrite  O' 

★ 

as  Q  in  (3.7)  and  (3.8).  Finally  we  can  get  the  same  representations  of  the 
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log-likelihood  ratio  test  statistic  as  given  in  (3.18),  (3.23) , (6. 15)  and 

(6.43).  By  the  same  way,  we  can  prove  the  analogues  to  Theorem  4.1,  6.1  and 

6.2. 
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